Backlund flux-quantization in a model of 
electrodiffusion based on Painleve II 



A.J. Bracken 1 , *L. Bass 1 Wd C. Rogers 2,3 * 
1 Department of Mathematics 
The University of Queensland 
Brisbane, Australia 
2 Department of Applied Mathematics 
The Hong Kong Polytechnic University 
Hong Kong 

3 Australian Research Council Centre of Excellence 
for Mathematics & Statistics of Complex Systems 
School of Mathematics and Statistics 
The University of New South Wales 
Sydney, Australia 



PACS Nos: 02.30.Hq 02.20.Bb 82.39.Wj 87.16.D- 87.10.Ed 

Key words: Backlund transformations; Painleve II equation; flux-quantization; 
electrodiffusion; liquid junctions; Nernst-Planck system; nerve conduction 

Abstract 

A previously-established model of steady one-dimensional two-ion 
electrodiffusion across a liquid junction is reconsidered. It involves 
three coupled first-order nonlinear ordinary differential equations, and 
has the second-order Painleve II equation at its core. Solutions are 
now grouped by Backlund transformations into infinite sequences, par- 
tially labelled by two Backlund invariants. Each sequence is charac- 
terized by evenly-spaced quantized fluxes of the two ionic species, and 
hence evenly-spaced quantization of the electric current-density. Fi- 
nite subsequences of exact solutions are identified, with positive ionic 
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concentrations and quantized fluxes, starting from a solution with 
zero electric field found by Planck, and suggesting an interpretation 
as a ground state plus excited states of the system. Positivity of 
ionic concentrations is established whenever Planck's charge-neutral 
boundary-conditions apply. Exact solutions are obtained for the elec- 
tric field and ionic concentrations in well-stirred reservoirs outside 
each face of the junction, enabling the formulation of more realistic 
boundary-conditions. In an approximate form, these lead to radia- 
tion boundary conditions for Painleve II. Illustrative numerical solu- 
tions are presented, and the problem of establishing compatibility of 
boundary conditions with the structure of flux-quantizing sequences 
is discussed. 
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1 Introduction 



Transport of charged ions across liquid junctions plays a fundamental role in 
a variety of important physical and biological contexts, in particular semi- 
conductor theory, electrochemistry and models of the nervous system. Many 
years ago, an extension of the Nernst-Planck model of ion transport [U, |2] 
was proposed IH E], incorporating the effect of the electric field that de- 
velops within the junction in response to charge separation. Such separation 
tends to develop as ions diffuse with different diffusion coefficients, while the 
consequent electric field acts to inhibit it. This feedback leads to nonlinearity 
of the extended model. 

Mathematically, the extension involves the use of Gauss' Law to relate 
charge density to electric field, together with the imposition of Einstein's 
relation between diffusivities and mobilities of the ions. In its simplest form, 
the model deals with one-dimensional transport across an infinite slab occu- 
pying < x < 5, for two types of ions carrying equal and opposite charges. 
With the concentrations of the two ionic species and the induced electric 
field within the slab denoted by c±(x) and E(x) respectively, the governing 
system of coupled ordinary differential equations (ODEs) obtained is 

c+ \x) = {ze/k B T) E{x) c+(x) - <$>+/D + , 
c_'(x) = -(ze/k B T)E(x) c_(x) - , 

E'(x) = {Anze/e) [c+(x) - c_(x)] (1) 

for < x < 5. Here $ + and $_ denote the steady (constant) fluxes of 
the two species in the x-direction across the slab, z their common valence, 
and D± the corresponding diffusion constants, while k B denotes Boltzmann's 
constant, e the electronic charge, and T the ambient absolute temperature 
within the solution in the slab. An important auxiliary quantity is the electric 
current-density 

J = 5e($+ - $_). (2) 

The slab equations ([I]) when supplemented by two-point boundary con- 
ditions (BCs) specify the electrical structure of liquid junctions including 
nerve membranes. The latter are so thin (of the order of 10~ 6 cm) that typ- 
ical physiological potential differences of 0.1 V generate large fields of order 
10 5 V/cm. In that context it is not only the potential difference across the 
slab (which was of primary interest to Nernst [T] and Planck [2]) but the 
electric field distribution within it that may be important [6]. 
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In a special case when $_ = 0, it was shown [3] that as a consequence of 
(II]), c_(a;) satisfies a second order nonlinear equation that can be transformed 
into the Painleve II ODE (PII). Independently, it was shown [3] that in 
the general case, and with 2-point BCs, the field E(x) also satisfies a form 
of PII (see (13) below). These were perhaps the first applications of PII 
to physics or biology. Since the introduction of the model defined by (jl|, 
this aspect and others have been much studied (see [7J IH |9] and references 
therein). In particular, more recent analysis [TQJ, [TT] has drawn attention to 
the relevance to the model of already-known exact solutions of PII, and of 
Backlund transformations p2] of its solutions more generally. 

In the formulation given in j3], the model equations (jl]) were supplemented 
by BCs corresponding to charge-neutrality at each face of the slab, so that 
c+(0) = c_(0) and c+(S) = c_(5). From these imply that E'(0) = = 
E'(5), and hence lead to BCs of Neumann type for PII as it appears in (13) 
below). We shall discuss the relevance of these BCs, and the possibility of 
imposing others corresponding to the presence of well-stirred ionic reservoir 
solutions outside the membrane slab, occupying — oo < x < and 5 < x < 
oo. We shall also downplay the role of PII, emphasized in previous analyses 
of the model, in part because of the peculiar feature that boundary values 
of E(x) appear in the coefficients of that ODE for E(x) itself, as it appears 
in the form (13). This complicates its theoretical and numerical analysis 



In contrast, we find that the system ([I]) is solved numerically with 
no difficulty by commercial packages such as MAPLE |15j . at least in cases 
corresponding to charge-neutrality at x = and x — 5. (See Sec. 6.1.) 

By working directly with the system ([!]), we also avoid singling out the 
electric field, and give comparable weight to the ionic concentrations, and 
their associated fluxes. This will enable us to elucidate the role of Backlund 
transformations in producing sequences of solutions of ([I]) in which the fluxes 
are quantized. 

Of particular interest are finite subsequences of exact solutions with pos- 
itive concentrations, generated by Backlund transformations from a solution 
of 0, first written down by Planck [2], with c + (x) = c_(x) and E(x) = 
throughout the junction. It is remarkable that Planck's solution, obtained 
by him when modelling electrodiffusion 10 years before his fundamental work 
on quantization of black-body radiation, should now appear some 120 years 
later as the seed solution in sequences exhibiting Backlund flux-quantization. 

As we shall see, there are difficulties reconciling either of the two kinds 
of BCs mentioned above, with the formation of such sequences of solutions. 
This raises the problem of identifying what form of BCs if any is consistent 
with 'Backlund flux-quantization.' 
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2 Preliminary remarks and Painleve II 

With the introduction of a suitable constant reference concentration c re f_, to 
be identified later, and with the change to dimensionless independent and 
dependent variables 

x* = x/5 , c±(x*) = c±(x)/c ref . , 

E*(x*) = (zed/k B T)E(x) } (3) 
and dimensionless constants 



At = -$ ± 5/c ref D ± , A = ^/ek B T/4ir(ze)W c ref . (4) 
the equations Q become (on immediately dropping all the asterisks) 

c + '(x) = E(x) c+(x) + A + , 

c_ '(x) = -E(x) c-(x) + A_ , 

X 2 E'(x) = c + (x) - c-(x) (5) 

for < x < 1. The current density is also conveniently made dimensionless 
by setting 

3= (~ze)c ref XD + + D_) J - (6) 
In terms entirely of dimensionless quantities we then have 

j = a-A_ - a + A + , (7) 

where 

a± = D±/(D + + DJ) (so that a± > and a + + a_ = 1) . (8) 

After the first two of @ are added, and the third is applied, one integra- 
tion can be performed to deduce that 

P{x) — 6x = B (const.) , (9) 

where 

P(x) = c+(x) + c_ (x) - \\ 2 E(xf , 9 = A + + A_. (10) 



Eq. ([9]) has the general meaning of a work-energy relation. This can be 
seen by returning to the original variables for the moment and multiplying 
through by k B T, to obtain the result 

k B T[c + (x) + c_(x)] - eE(x) 2 /%n - + = const. (11) 

where u + = D + /kBT and ii_ = D_/k B T are the mobilities of the two types of 
10ns. Here k B T[c+(x) + c_(x)] - eE(x) 2 /8 tt is the ionic (osmotic) pressure at 
x in the presence of the electric field. The minus sign reflects the fact that the 
electric field tends to keep the oppositely charged ions together, effectively 
reducing their degrees of freedom and hence the associated pressure. The 
term [§ + /u + ]x (resp. [$_/tt_]x) is the work done per unit volume against 
the resistance offered by the background fluid, in transporting positive (resp. 
negative) ions from to x. Accordingly, we may refer to P(x) and 9 as (the 
dimensionless forms of) the pressure at x, and the resistive force-density, 
respectively. 

Differentiating the third of ^ once and then applying the first two of 
^ to eliminate c + '(x) and c_ '(x) from the result, leads to 

\ 2 E"{x) = E(x) [c + (x) + c_(z)] + (A + - A_) (12) 

and, with the help of (J9| to 

X 2 E "(x) - \\ 2 E(xf - (6x + B)E(x) - (A+ - A_) = . (13) 

This can be brought to a standard form of PII by a constant scaling of the 
dependent variable, together with a linear transformation of the independent 
variable [HE]. 



From any solution E(x) of (13), corresponding c±(x) can be obtained 
using g and the third of @. In this way, PII determines all solutions of 
([5]). The results of many studies [ID]- [II], [IS]- [21] of PII can then be brought 
to bear on the problem of interest, in particular Backlund transformations 
and exact solutions. These have been extensively explored elsewhere [TT] in 
the context of electrodiffusion on a half line rather than a finite interval of 
the x-axis, which is our focus here. 

3 Planck's solution 

Planck [2] considered the situation where A is vanishingly small. As Q shows, 
this occurs when the electric charges ±ze are much larger in magnitude 
than ^ekgT /And 2 c re f.. (Here the choice of c re j. becomes relevant — see the 
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Remark following (67) below.) In the limit A — > 0, with c(0) = Co, the system 



^ is easily solved to give Planck's solution 

c + (x) = c_(a;) = c(x) = c + \{A + + A_)x , 

E(x) = l(A_-A + )/c(x). (14) 
With the help of @ the electric field may also be written as [5] 

E(x) = ^- + {a+ - a - ) ^- Co) , (15) 

where c\ = c(l) = cq + (A + + A_)/2, which brings to the fore the role of 



the electric current density: the first term on the RHS of (15) is the ohmic 
contribution to the field, and the second term is the gradient of the Nernst- 
Planck potential [2J. The situation where j = is of particular physical 
interest. 



For any value of j, Planck's solution (14) might be used as the zeroth- 
order term in an asymptotic expansion, in powers of A 2 , of a solution to 
the system ^ with A / 0. This would require the methods of singular 
perturbation theory, as is clear from the way in which A 2 appears multiplying 



the highest derivative in (13). That approach has been developed previously 



for closely related problems [251 126] . and will not be pursued here. 

For our purposes in what follows, it is more important to note that ( 14 ) 
defines an exact solution of the full system ^ in the special case that A + = 
A_ = ci — c and hence E(x) = 0. This solution, which has j ^ 0, plays the 
role of the seed for an infinite sequence of exact rational solutions [HJ UHl EH 
|22| [TUl ITT] generated by Backlund transformations, as highlighted in Sec. 6 
below. 



4 A discrete symmetry group 

The system ^ must be supplemented by appropriate BCs, as we discuss 
below, to determine physically meaningful solutions with nonsingular c±(x), 
positive on [0, 1]. For the time being we ignore such considerations. Bearing 
in mind the way that PII determines all solutions of ([5]), and that solutions 
of PII itself, when regarded as analytic functions, may have (movable) sim- 
ple poles, we mean by "a solution S = (c+ , c_ ,E , A + , A_)" , any two real 
constants A + , A_, together with any three real- valued functions c+, c_, E 
satisfying (fsl) almost everywhere on [0,1]. 
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Given one solution S = (c+ , c_ , E , A + , AJ), straightforward if some- 
times lengthy manipulations confirm that in general there exist four other 
solutions. These transformed solutions can be found from known symmetries 
of PII [H], but they are presented together here, in terms of the variables of 
the system ([5]), for the first time. The transformed solutions, each of which 
we denote by S = (c+ ,d-,E, A + , A J) in turn, are as follows :— 

4.1 Conjugate solution 

Here S = C(S), where 

c+ = c_ , c_ = c+ , E = -E , A+ = A- , A- = A + . (16) 

Noting that a second application of C returns S to the starting solution 
S, we write 

C 2 = I. (17) 

4.2 Reflected solution 

Here <S = H(S), where 

c+(x) — c + (l — x) , c_(x) = c_(l — x) , 

E(x) = -E(l - x) , A+ = -A+, A_ = -A_. (18) 
We note that 

k 2 = i, nc = cn. (19) 

4.3 Backlund-transformed solution 

Here S = B(S), where 

c+ = c_ + 2\ 2 A + E/c + + 2\ 2 A 2 + /c 2 + , c_ = c+ , 

E = -E- 2A + /c + , A+ = 2A + + A_, A_ = -A + . (20) 
We note that if A + = 0, then B(S) = C(S). 
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4.4 Inverse Backlund-transformed solution 

Here S = B~ 1 (S), where 

c+ = c_ , c_ = c+ - 2\ 2 A_E/c_ + 2\ 2 A 2 _ /c 2 



E = -E + 2A-/c- 



A, 



-A- , A_ = 2A_ + A, 



(21) 



We note that if A_ = 0, then ^(S) = C(5). 

It is straightforward to check that B~ l is indeed inverse to B, and we 
write 



-i 



B~ L B = I = BB 
It is also easily checked that 

cb = b- 1 c, nB = Bn, izb- 1 = b- x iz 



(22) 



(23) 



From (23), (19) and (17), it then follows that C, 1Z and B, together with 



their inverses and with I as unit, generate a discrete invariance group of 
transformations of the system Except on special solutions, this group is 
of infinite order. 



From (16), ( |20| and (21), it is also easily seen that the quantities P(x) 
(pressure at x) and 9 (resistive force density) as in @ and p| are invariant 
under the action of the subgroup generated by C, B and B~ x . It is unclear 
why these physical quantities should be singled out to be invariants. The 
invariant 9 plays a central role in flux-quantization below. 



4.5 Gambier-transformed solutions 

There also exist invertible transformations, which we call Gambier and in- 
verse Gambier transformations, that carry solutions back and forth between 
two special subclasses. These subclasses are of particular interest because 
one contains, among others, exact solutions of rational type, while the other 
contains, among others, exact solutions of Airy type These transforma- 
tions and their inverses, to be denoted by Q± and Q^ 1 ^ respectively, and to be 
defined in what follows, can also be derived from a known (restricted) sym- 
metry of PII discovered by Gambier [T?] |2"3"1 12"1] . We present them explicitly 
here in terms of the variables in ^ for the first time, as follows. 

The transformation Q + (resp. QJ) converts any solution S in which A + = 
(resp. A_ = 0), into a solution S in which A + = A_. If S has the value B 
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for the first-integral in ([9]), then <S has corresponding value B = —B/2. In 
detail, if A + = and A_ = A, the transformed solution S = G+(S) has 

c+(x) = \XE(x)y/2c+(x) + \c+(x) - \{Ax + B), 

c_(z) = -fA£(x)^2c + (x) + \c + {x) - \(Ax + B), 

E{x) = ^2c + (x) /A, A+ = i_ = -|A, B = -\B. (24) 
Similarly, if A_ = and A + = A, the transformed solution Q-{S) has 
c+(x) = -\\E(x)y/2c-{x) + |c_(x) - + 



C_ X 



iAE(x)v/2c_(x) + |c_(x) - 1(Ae + B) 



= y/2c-(x) /A, i+ = i_ = -|A, B = -\B. (25) 

4.6 Inverse Gambier-transformed solutions 

These convert any solution S for which A + = A_ = A, with a given value of 
B, into solutions with A + = (resp. A_ = 0) and with B = —2B. In detail, 
S = G+\S) has 

c+(x) = \\ 2 E(x) 2 , £(x) = 2(c+(x) - c_(x))/A 2 £(x) , 

c_(x) = -\\ 2 E(xf + 2(c + (x) - c4x)) 2 /X 2 E(x) 2 - AAx - 2B , 

i+ = 0, i_ = -4A, B = -2B. (26) 
and similarly, ^Z 1 (5) has 

c_(x) = \\ 2 E(xf , = -2(c+(z) - c_(x))/A 2 E(x) , 

c + (x) = -^A 2 £(x) 2 + 2(c+(x) - c_(x)) 2 /A 2 £(x) 2 - 4Ar - 25 , 

i_ = 0, i+ = -4A, B = -2B. (27) 

The case when i?(x) = needs special treatment (see Appendix A). Note 
that 

GZ\S)=CG-\S). (28) 
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5 Sequences of solutions 

Repeated application of B and B~ x produces from any 'seed' solution S, a 
sequence Qs of solutions that, except in special cases (see below), is doubly 
infinite. Thus 

Qs = {..., B-\S) , B-\S) , S , B(S) ,B 2 (S),...} (29) 
There is an associated conjugate sequence 

Q c{s) = {. . • , B- 2 C(S) , B- X C{S) , C(«S) , BC(S) , B 2 C(S) , . . . } 

= C{... ,B 2 {S),B{S), S,B-\S), B- 2 (S),...} 

= CQ T S , (30) 

where Q$ is the transpose of the sequence Qs, i.e. the sequence obtained 
from Qs by replacing the nth member by the (— n)th member, for n = 
, ±1 , ±2 , . . . . There is also an associated reflected sequence 

Qns) = {... ,B- 2 n(S),B- 1 n(S),n(S),Bn(S),B 2 n(S),...} 

= K{..., B- 2 (S) , B~\S) , S , B(S) , B\S) , . . . } 

= HQs. (31) 

Because the quantities 9 and P(x)—6x (with value B) are invariant under 
the action of C, B and B~ l , every solution in Qs and Qc(s) has the same 
values for these quantities. On the other hand, in the reflected sequence 
Qn(S), every solution has B = B + 6 and 6 = —6. In this way, each of the 
sequences Qs, Qc{s) an d Qn(s) is partially labelled by the values of B and 6. 

6 Backhand flux- quantization 



It follows from (20) and ( |2l[ ) that, if we write for the moment 

B n (S) = (cJ\cL n \£ (n \Af\AL n) ), n E {0 , ±1 , ±2 ,...}, (32) 

then 

A (n) + ^(n) = 9 
A (n) _ A (n) = A (0) _ A (0) + ^ (33) 
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Here 9 may take any real value. If 9 ^ 0, we see from (33) that without loss 
of generality we can assume for S that 



and we then have 

A (n) 



A 



(0) 



.4 



(0) 



-9 < (f) < 9 



(34) 



a + n)9 + U 



A 



in) 



n)9-U, n E {0 , ±1 , ±2 ,...}. 



(35) 



It follows that the steady ionic fluxes as in Q are also quantized in each of 
the sequences of solutions generated by Backhand transformations and their 
inverses, with values for the n-th element of the sequence given by 



4> 



(n) 



n + 1)$^ 0) +n (£>+/£>_)$ 



(0) 



4> 



(n) 



(n - 1)$ 



(0) 



n(D_/D + )<f> {0) 



+ > 



(36) 



and then the electric current density as in (|2| is also quantized, with n-th 
value given in dimensionless form by 



J 



(n) 



- a+(0 + 4>)]/2 - n9(a + + a_) 



and in dimensional form by 



J{ n) 



ze 



[n + l + n(D_/D + )}^ ] + [n - 1 + n(D + /D_)]^ 



J {0) +nAJ, A J = ze{D + + _D_ 



(0) 



(0) 



(37) 



(0) 



(3f 



The case 9 = (zero resistive force-density) is degenerate; then may take 



any real value, and there is no flux-quantization, as (35) shows. In this case 



the system (|5j) is solvable explicitly in terms of elliptic functions, correspond- 
ing to the reduction of PII in Q to an elliptic ODE [21 HH E3- 

Flux quantization is a remarkable feature of the structure of solutions of 
((!]) when grouped into sequences by Backlund transformations. To begin to 
explore the physical interpretation of this mathematical result, we reconsider 



Planck's electrically neutral solution (14) in the case E(x) = when it 
becomes an exact solution of (|5]), that is, when A + = A_ = c\ — Co = A. 
We take this as seed solution S in a doubly-infinite sequence Qs of exact 
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solutions generated by Backlund transformations and their inverses. Writing 
for the moment <S = B(S) and = B~ 1 {S) for the two members of this 
sequence on either side of S, we find using pOl and (f2lT) that S has the form 



c + (x) = c + Ax + 2X 2 A 2 /(c + Ax) 2 , c_(x) = c + Ax . 



E{x) = -2A/(c + Ax) , A + = 3A, A_ = -A; 



(39) 



S has the form (from (14)) 

c + (x) = c_(x) = c + Ax 



E(x) = , A + = A^ = A; (40) 

and «S' has the form 

d(x) = c + Ax + 2X 2 A 2 /(c + Ax) 2 , c^(x) = c + Ax , 



E\x) = 2A/(c + Ax) , A\ = -A, A ] = 3A 



(41) 



We note that (|4T| can also be obtained from (39) using (23). In each of 



these three solutions, the concentrations are everywhere positive as required. 
It is natural to think of S as describing a 'ground state' of the system, where 
E(x) = and j = (ct_ — a + )A. Then <S and describe 'excited states,' 
with j = — (3a + + a~)A and j = (a + + 3a^)A respectively. Note however 
that charge-neutrality at x = and x = 1 is satisfied by the seed solution S, 
but not by S and cv, and it is not clear what interpretation to give to the 
boundary-values satisfied in these two excited states. 

The sequence Q$ generated from S in this case corresponds to the well- 
known doubly-infinite sequence of rational solutions of PII [321 ESI EB E21 
ITUl ITT] . Expressions for the members of the sequence beyond those given in 



(39)-(41) become increasingly more complicated, but it is straightforward to 



use MATLAB [28J to construct such members approximately without solving 
(5| numerically, simply by proceeding from the seed solution using (20) and 
@). 

Fig. 1 shows c + (x) as in (40) and its first eight iterates under B, for the 
choices Co = 1/3, A = 1/3, 



and A 2 = 0.01. Note from Q that these are 
also the graphs of the first nine iterates under B of c_(x), and from (21), 
they are also the seed and first eight iterates under B~ l of c_(x), and also 
the first nine iterates under B~ l of c + (x). The ninth graph shows unphysical 
negative concentration values, showing that for these parameter values, there 
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Figure 1: On the left, graphs of c+(x) (straight, solid) as in (40) and its nth 
transforms under B, for n — 1 (solid), n = 2 (dashed), n = 3 (dotted) and 
n = 4 (dashdotted). On the right, graphs for n = 5 (dashdotted), n = 6 
(dashed), n = 7 (dotted) and n = 8 (solid), showing appearance of unphysical 
negative values in the final case. 



is a ground state plus seven excited states plus their seven conjugates. The 
corresponding evenly-spaced current-density values, from the seventh iterate 
under B~ l through to the seventh iterate under B are, in dimensionless and 
dimensional form 

j (n) = -[(2n + l)a+ + (2n-l)a_]A, 



= ze(c - Cl )[(D+ - D_) + 2n{D + + D_)]/5, 



n = -7, -6, •••,7. (42) 

Experiment with this technique using MATLAB shows that the number of 
excited states with positive concentrations obtained from (14) (in the case 
with E(x) = 0), depends on the values of c , C\ and A. 

Examples like this show that for any physically sensible seed solution 
S, we must expect that only a finite subsequence of Qs may correspond 
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to physically meaningful situations; in other words, only a finite number of 
physically acceptable solutions will be generated from any acceptable seed 
solution by Backlund transformations and their inverses. Within that subset 
of solutions, however, we have nevertheless the remarkable result that ionic 
fluxes and the electric current-density are quantized in accordance with (36). 

A notable feature of (38) and (42) is that increments of fluxes pertaining 
to successive values of n are all equal and all determined at n = 0, despite 
radical differences between the solutions of ([!]) associated with those fluxes 
at different n values. What common element of such disparate solutions can 
engender consistency with the equality of the flux increments? We develop 
elsewhere [29] the hypothesis that the common element is the quantum of 
electric charge. 

[Remark: In regard to the general question as to when Backlund or inverse 
Backlund transformations of a given physically meaningful solution of ^ 
give rise to new (potentially) physically meaningful solutions, with positive 
concentrations, we may make the following intuitive observations. When 
D + ^ it is clear that no steady electric field can counteract diffusion of 
both ionic species and force both fluxes $ + and $_ to zero. However, there 
is a field capable of equating the fluxes, so that positive and negative ions 
moving together make up zero electric current density; this is the situation 
where J — in ([2]). To achieve this, E(x) must be directed so as to oppose 
the flux of the more mobile species. Because E(x) has to be anti-parallel 
(resp. parallel) to the flux of positive (resp. negative) ions in order to 
oppose it, if D + > D_ then E(x)& + < 0, and if D + < D_ then E(x)&_ > 0. 
In dimensionless form, if j = and a + > «_ then E(x)A + > 0, and if 
j = and a + < a_ then E(x)A- < 0. But then it is clear from (20) and 
( |2~lj ) that if c + (x) and c_(x) are positive, so are their Backlund transforms 
c+{x) and c_(x) in the first case, while their inverse Backlund transforms 
are positive in the second case. Note that this result is independent of the 
choice of boundary conditions. For j ^ the argument breaks down; and 
since a solution with j = is carried by a Backlund or inverse Backlund 
transformation into a solution with j '• ^ in general, as (37) shows, the 
argument cannot be applied more than once to a solution with j = 0; then 
we cannot easily tell if further Backlund or inverse Backlund transformations 
will produce solutions with positive concentrations.] 

It is clear that an analysis of BCs and positivity requirements for solutions 
of ^ is called for; we turn to this in the next section. 
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7 Boundary conditions and positivity of con- 
centrations 

The ionic fluxes $+ and appearing in 0, and hence the dimensionless 
constants A + and A_ appearing in d5|, are not to be considered in general as 
given. The electric current density Q and hence the linear combination ^ 
may be prescribed, but this leaves at least one of A + and A_ to be determined 
as part of any solution S = (c+ , c_ , E , A + , A_) of dq). Five pieces of data 
must be supplied to fix such a solution. 

So far we have partially characterized the seed solution S of a sequence 
Qs by the values B, 9 and 0, thus fixing the values of A + and A_. Ignoring 
for the moment that this may not be possible in practice, we note that the 
conjugate sequence Qc(s) is then also partially characterized, by correspond- 
ing values B, 9 and — 0, and the reflected sequence Qn(S) is also partially 
characterized, by corresponding values B — 9, —9 and — 0. Two further num- 
bers are needed to fix each of these three sequences uniquely. They can be 
taken to be the two further numbers that are needed to fix the seed solution 



S = (c + , c_ , E , A + , A_), which so far has, using (J9J) and (35), 
A + = (9 + <P)/2, A_ = (0-0)/2, 

c+(0) + c_(0) - ±\ 2 E(Q) 2 = B = c+(l) + c_(l) - ^A 2 £(l) 2 - . (43) 

For example, we could prescribe the values at x = of all three of c+, c_ 
and E, bearing in mind that we have already prescribed the value B of the 
combination c+(0) + c_(0) - X 2 E(0) 2 /2. Then we can reasonably hope that 
the seed solution S is uniquely determined, and so in turn, every element of 
the sequences Q s , Qc(s) and Qn(s)- 

Writing again B n (S) = (c+ , c_ , E , A + , A_), we could then determine 



the values of c+(0), c_(0) and £"(0) by repeated use of the formulas (20) 



and (21) at x = 0, although the resulting expressions quickly become very 
complicated as \n\ increases. Similarly, if we were to fix S by prescribing the 
values of c + (l), c_(l) and £7(1), we could then determine the values of c + (l), 



c_(l) and E(l) by repeated use of the formulas (20) and (21) at x = 1. Note 



however that if we were to fix the seed solution by giving the values of, say, 



E(0) and E(l) in addition to M3 



then we could not determine the values 



of E(Q) and E(l) using (|20j) and (|2lJ) alone. 

Because the values of A + and A_ are not in general known a priori, such 
considerations are of more mathematical than physical interest, and we now 
turn to more physically relevant ways of prescribing the five pieces of data 
needed to fix a solution of Two types of BCs are of particular interest. 



14 



7.1 Charge-neutral boundary conditions 



These were adopted by Planck [2] and in the formulation of the present model 
given in [U E], and take the form 



c+(0) = c_(0) = c > and c+(l) 



;i) = Cl >o. 



(44) 



Without loss of generality, we can suppose that cq < c\. Taking c re /. to equal 
the sum of the concentrations at the two faces of the slab before proceeding 
to dimensionless variables as in ([3]), we then have in dimensionless form 



< c < Ci , c + Ci = 1 . 



(45) 



Eqs. (44) provide four of the five pieces of data that we expect are needed 



to determine a solution of The remaining piece is typically provided by 
fixing the value of the current density, 



3 = Jo ■ 

Using the last of ([5j , we see from ( 44 ) that 

E'(0) = = E'(l) 



(46) 



(47) 



implying Neumann BCs for PII as in (13). However this boundary- value 



problem for ( 13 ) is not of a standard form because constants to be determined 



as part of the solution appear in the coefficients of the ODE itself. In full, 
we have using ([9]), (44) and (46) in (13) that 



X 2 E"[ 



x 



\\ 2 E{xf 



[2c - \\ 2 E(0) 2 + {2( Cl - c ) + l\ 2 (E(0) 2 - E(l) 2 )} x] E( 



x 



— (a+ — a- 



{2( Cl - c ) + l\ 2 (E(0) 2 - E(l) 2 )} - 2j • (48) 



Here A, Co, ci, a+, a_ and jo are known constants, while E(0 and E(l) are to 
be determined together with E(x), for < x < 1. Despite its unusual form, in 
the case j = existence of positive, strictly decreasing solutions E(x) of (48) 
and (47) satisfying coE(0) < ciE(l), and non-existence of negative solutions, 
has been shown [13l HI] for a + > a_. By charge- conjugation, existence of 
negative, strictly increasing solutions satisfying cqE(0) > C\E(1), and non- 
existence of positive solutions, has been shown for a + < a„. If a + = a_, 
then jo — implies that A + = A_ = 0, and the explicit solution E(x) = 



exists as in (14) 
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From any of these solutions E(x), a solution to (JHl), (44) and ( [46] ) can be 
constructed using ^ and the last of Note that when E{x) is strictly 
increasing (resp. decreasing), then c+(x) > c_(x) (resp. c+(o;) < c_(x)) on 
(0,1). 

Uniqueness has not been established for any of these solutions, and there 
are no general existence or uniqueness results when jo ^ 0. 



It has proved difficult [14J to solve (48) directly by numerical methods 
in order to obtain approximate solutions that illustrate the results obtained 
analytically [T31 E] . This difficulty, which stems from the fact that the ODE 
to be solved involves constants that have to be determined as part of the 
solution, can be circumvented simply, by working with the system ^ rather 
than its consequence (48), or more precisely, by rewriting ^ as a system of 
five coupled first-order ODEs, 

c + '(x) = E(x) c+(x) + A+(x) , c_ '(x) = —E(x) c-(x) + A^(x) , 

\ 2 E\x) = c+(x) -c-(x), A+'(x)=0, A_'(x) = 0, (49) 
which does not involve any unknown constants. The MAPLE [15J routine 



dsolve solves (49) almost instantaneously when supplemented by BCs of the 



form (44) and (46), with the latter rewritten as 

a_A_(0)-a+A + (0)=j o . (50) 

Fig. 2 shows plots obtained in this way, for c = 1/3, c\ = 2/3, A = 0.5 
and jo = 0, showing expected behaviour [131 E] : a positive, monotonically 
decreasing E(x) on the left, in the case a + = 0.8 > a_ = 0.2, and a negative, 
monotonically increasing E(x) on the right in the case a + = 0.4 < a_ = 0.6. 
On the left, c E(0) w 0.15 < c x E(l) w 0.25; on the right, c E(0) w -0.049 > 
c x E{l) w -0.083. In each case, @ is satisfied. 

Comparison with the behaviour of E(x) in cases with j ^ 0, as in Fig. 3 
below, suggests strongly that these known existence and monotonicity results 
[131 114] can be extended to such cases also. 

It is easy to show as follows that for any solution of (|5]), (44) in which 



c + (x) and c_(x) are continuously once-differentiable on (0, 1), and continuous 
from the right (resp. from the left) at x = (resp. x = 1), then both c + and 
c_ are everywhere positive on [0, 1]. 

Suppose firstly that c + (x) < on some open subinterval of [0, 1]. Because 
c+(0) = Co > 0, and c + (l) = C\ > 0, it then follows that there is at least one 
point x\ and one point x<i in [0, 1] with 

c+(xi) = 0, c+'OO^O, c+(x 2 ) = 0, c + '(x 2 )>0. (51) 
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Figure 2: Graphs of E(x) when Co = 1/3, c\ = 2/3, jo = and A = 0.5. On 
the left, a + = 0.8, a_ = 0.2; on the right, a + = 0.4, a_ = 0.6. 

These conditions are only compatible with the first of (|5| if 

c + '(x 1 ) = c + '(x 2 ) = A + = 0. (52) 

But if A + — 0, it follows from the first of (jHJ and (|44f that 

c+ (x) = c e'° E & * , < x < 1 , (53) 

which is everywhere positive, hence providing a contradiction. 

Next suppose instead that c + (xs) = at some point x% G [0, 1]. Because 
it cannot happen that c + (x) < on an open subinterval of [0, 1], it follows 
that £3 is a minimum, and c + '(x^) = 0, so again from the first of ^ we have 



A + = 0, hence (53), and again a contradiction. 

Thus c + > on [0, 1], and similarly c_ > on [0, 1]. 

With that established, it is possible to prove for any solution of this type 
that at least one of A + and A_ must be positive if C\ > c . Assume to the 
contrary that A + < and A_ < 0. From ^ it then follows that 

d (e^ x) c+(x)) jdx = A + e^ {x) < , 
d (e-^ x) c_ (a;)) jdx = A. e~^ x) < , (54) 
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Figure 3: Graphs of c+(x) (dash), c_(x) (dashdot), E(x) (solid), A + (dot) 
and A_ (spacedot), corresponding to charge-neutral BCs pi] ), (46), with 
c = 1/3, Ci = 2/3, A = 0.7, a + = 0.4, a_ = 0.6 in each case, and with 
jo = 0.16 on the left, and jo = 0.6 on the right. Both A + and A_ are positive 
on the left; on the right, A_ is positive but A + is negative. 



where 

ij;(x) = - I E(y)dy. (55) 



Then 

e^°) c+(0) > e^^ 1 ) c + (l) e^ )"^ 1 ) > Cl /c 



-V(o) 



c_(0) > e"^ (1) c_(l) ^ e^ 0) -^ (1) < co/ci , (56) 



which are inconsistent with c\ > cq. 

Numerical experiments indicate that there are cases for which both of 
A + , A_ are positive, and cases where only one is positive. Fig. 3, obtained 
in the same way as Fig. 2 above, shows examples of both behaviours. Note 
that if jo — in (46), then rtTfi implies that A + and A_ have the same sign, 



and so both must be positive if c\ > cq. 
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We now show that BCs of the form (44) and (45) are not preserved under 
Backlund transformations. We have already seen this in the special case of 
the exact Planck solution, described in (14). Consider <S = B(S) where S 



is a solution satisfying (44) for some Co and c\ as in (45). From (20) we see 
that c+(0) = c_(0) requires that 



2A + \ 2 (c E(0) + A + )/cl = 



(57) 



or c+'(0) (= c E(0) + A+) = 0. Sim- 
0. But since S should be fixed by the 
five pieces of data in (|44|) and (|46|), it is not possible to impose any of these 



which in turn requires either 
ilarly, either A + = or c+'(l) 



conditions, and none will hold in general. Certainly there can be solutions S 



for which (44) hold together with A + = rather than (46), but then the cor- 



responding current-density j is not known a priori, but is determined as part 
of the solution. In such a case, B(S) reduces to C(S), the charge-conjugate 
of S, with A + = A_ and A_ = 0, as noted after (20). It follows that a second 
iteration of B is not possible unless it is also true that A_ = 0; but when 
A + = A_ = 0, no flux-quantization occurs, as (36) shows. Corresponding 
remarks apply for inverse Backlund transformations. Thus charge-neutral 



BCs (44) are not preserved in form under flux-quantization by Backlund 



transformations and their inverses. 



7.2 Corrected boundary conditions 



In proposing charge-neutral BCs (44), we have in mind reservoirs of well- 



stirred ionic solutions occupying — oo < x < and 1 < x < oo, with 
c + (x) and c_(x) both equal to Cq throughout the left-hand reservoir, and 
both equal to C\ throughout the right-hand reservoir, and with E(x) = 
throughout both reservoirs. Stirring leads to the diffusion coefficients D + 
and -D_ becoming effectively infinite, so that we also have A + = A^ = j = 
throughout both reservoirs (cf. Q and (J5|). When, as we expect, (44) and 



( |46) ) determine a solution of ([5| in < x < 1, boundary- values E(0 + ) and 
E(l_) of E(x) are determined in turn as part of that solution. There is no 
reason to expect that these values will be zero in general, and so they will not 
equal the assumed values in the reservoirs, -E(0_) = E(l + ) = 0. But such 
jumps in the values of E(x) across the faces of the slab imply the presence of 
nonzero surface charges, contradicting the assumption of charge neutrality 



there. From this point of view, the charge-neutral BCs (44) are inconsistent. 

In reality, we must expect that the ionic concentrations and the electric 
field are continuous across each face of the slab. Then, for the left-hand 
reservoir for example, c + (0) ^ c_(0) and E(0) ^ 0, but as x — > — oo, c+(x) 
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and C-(x) decay to a common positive value C(_ 00 ), while E(x) decays to 
zero. Similarly, in the right-hand reservoir as x — > oo, c + (x) and c_(se) decay 
to a common positive value C( +oc ), while E(x) decays to zero. Without loss 
of generality, we can suppose that C( +oc ) > C(_oo). It is now appropriate to 
choose the sum of these two limiting concentrations, before the change to 
dimensionless variables, as the c re f in (|3j) and Q; in dimensionless terms we 



then have in place of (45) that 



< c ( _ oo) < c (+00 ) , q.oo) + c (+00 ) = 1 . (58) 

In the left-hand reservoir, equations of the form ([T| apply in simplified 
form, as D + — > oo and D- — > oo; in dimensionless form, we get equations of 
the form ^ with A + = A- = 0, so that for — oo < x < we have to solve 

c + \x) = E(x)c + (x) , c_ \x) = —E(x)c~(x) , 



\ 2 E'(x) = c+(x) -c_(x), (59) 

with c±(x) — > C(__oo) and E(x) — > as x — )■ — oo. Then the first two equations 
integrate to 

c + (x) = c hoo) e-vW ; c _( x ) = C( _ oo) e ^) ; ( 60 ) 

where 

/X 

E(x) dx (61) 
-oo 

is the electrostatic potential of the electric field in the reservoir. Because 



E(x) = —(p'(x), the third of (59) then gives the nonlinear Poisson-Boltzmann 
equation 

- AV (a;) = c ( _ oo) [e-*W _ e <**)] , (62) 

to be solved on — oo < x < subject to ip(x) — > and v 9 '( x ) — ^ as 
x — >■ — oo. The exact solution of this problem is known [301 E]- In Appendix 
B we present this solution and describe the highly nonlinear exact BCs at 
the junction face(s) to which it leads. Here we suppose that 

\<p(x)\ «1, (63) 



so that ( 62 ) linearizes to 

^(f/\x)=(p(x), Ao = A/v^^. (64) 
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with solution matching the limiting BCs at x — — oo given by 

(p(x) = ip(0) e x/Xo . 



(65) 



In this approximation we then have, from (60) and (61) 

c 



_(x) = c ( _ oo) [l + <p(0)e x ^] , 



E{x) 



y(Q) 

An 



3 a:/Ao 



(66) 



showing that An (in dimensional form, An<5) can be interpreted as a Debye 
shielding length for the reservoir-junction interface. It is the characteristic 
length over which c+(x), c_(x) and E(x) decay to their limiting values as x 
decreases from the interface value x = 0. In terms of the original variables, 



X S = yj ek B T/8Tr(ze) 2 C(- 



oo) 



(67) 



[Remark: Just as Ao is a Debye length at concentration c_oo, so the con- 
stant X/y/2 of Q is a Debye length at concentration c re f.. In analyzing 
mathematically situations where A approaches zero [32] (see Sec. 3), a dif- 
ferent choice of c re f. in Q may be more appropriate than the one we have 
made at the start of this section. For example the choice of the logarithmic 
mean (c( +oc ) — C(_ 00 ))/ln(c( +00 )/c(_ 00 )) as c re /. has the advantage that c re /. 
then tends to zero like the smaller of C(_oo) and C( +00 ), guaranteeing that as 
A approaches zero, the Debye length in each reservoir also approaches zero.] 

In (66), the value of (p(0) is determined from the solution of ^ in the 
slab, which fixes c±(0) and E(0). 



From (66) we have 



c + (0) + c_(0) =2^), 



c+(0) - c_(0) = -2c ( _ oo) ^(0) = 2c ( _ oo) A £(0) 



and then, from the third of equations ( 59 ) 



\ E'(0) = E(0). 
Similar considerations in the right-hand reservoir lead to 



(6* 



(69) 



c + (l) + c_(l) = 2c 



(+oo) 



c+(l)-c_(l) = 2c (+oo) A 1 E(l) 



(70) 
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and 



\ 1 E'(1) = E(1) 



(71) 



In (70) and (71), Ai is defined like Ao in (64), with C( +oc ) replacing C(_oo). 

With C(_oo), C( +00 ) and j given constants, equations (68), (70) and (46) 
provide the new BCs that are to be applied to fix a solution of the system 
x < 1, replacing MM and (|46[) . Fig. 4 provides an illustration 
of the solution to the combined problem (left reservoir, slab, right reservoir). 
In the slab, equations ^ were solved numerically using MAPLE, with the 
BCs (pi) and ((701. 



Figure 4: Graphs of c + (x) (dash), c_(x) (dashdot), E(x) (solid), A + (dot) 
and A_ (spacedot), in the left reservoir, in the slab, and in the right reservoir, 
corresponding to corrected (approximate) BCs. Parameter values are C(_oo) = 
1/3, c (+oo) = 2/3, A = 0.7, a+ = 0.4, a_ = 0.6, with j = 0.16. 



This boundary value problem can also be posed in terms of PII, which in 



this case takes a simpler, more standard form than (48). To see this, note 
firstly that Co and c\ in (48) are now to be replaced by C(_oo) and C( +00 ), 



respectively, according to the first of (68) and the first of (70). Then note 



that, as a consequence of (63) and the third of (66) 



\\ 2 E(0) 2 < 2c ( _ 



(-0O) 



(72) 



Similarly 



^(lf^q 



+oo) 



(73) 
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The terms on the left-hand sides of these inequalities are therefore to be 
neglected relative to the terms on the right-hand sides in the ODE that 



replaces (48), and it therefore takes the form 



\ 2 E"(x) = l\ 2 E(x) 3 + 2[c { ^ oo) + {< 



'(+00) 



"(-do) 



}x] E( 



-2(a+ - a_) {c( +oo) - C(_oo)} -2j 



(74) 



This equation is to be solved on < x < 1, subject to the radiation BCs 
(69) and (71). The constants C(_oo), C( +oc ), a+, a_ and jo appearing here are 
to be regarded as given. Nothing is known about existence or uniqueness of 
solutions to this problem, to our knowledge, nor about positivity of resultant 
concentrations. 



Given (72) and (73), the question arises if \\ 2 E(x) 2 should also be ne- 
glected relative to C(_oo) and C( +oc ) for all < x < 1, and not only at the 



boundaries. If so, then the cubic term in (74) drops out, and the resultant 



linear boundary-value problem should be exactly solvable in terms of Airy 
functions (cf. [3]). This will not be pursued here. 

Similar manipulations to those in the preceding subsection show that, 



like charge-neutral BCs, those in (68), (70) and (46) are also not preserved 



in form under Backlund transformations and their inverses. 



8 Concluding remarks 

Quantization of ionic fluxes and of the electric current-density by Backlund 
transformations of solutions, in this one- dimensional model of electrodiffu- 
sion within a liquid junction, is a remarkable mathematical phenomenon. 
What physical significance can be attached to this phenomenon is not com- 
pletely clear, but the example provided by the exact Planck solution and 
its transforms is highly suggestive of a corresponding ground state and ex- 
cited states of the system. A difficulty with this interpretation, and with the 
concept of flux-quantization by Backlund transformations more generally, is 
that BCs that are invariant in form under such transformations have not yet 
been identified. 

We have seen that familiar charge-neutral BCs, first introduced in the 
context of electro-diffusion by Planck, lead to several attractive properties 
of solutions. Existence of such solutions has previously been proved under 
certain conditions [131 E], and we have shown positivity of the functions 
representing ionic concentrations. However, we have argued that such BCs 
are not physically consistent. Moreover, they are not preserved in form under 
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Backlund transformations. This has prompted us to examine more realistic 
BCs. To find these we have had to solve for the exact form of ionic concen- 
trations and the electric field in well-stirred reservoirs outside each face of the 
junction. The new BCs define novel boundary-value problems for the system 
of ODEs that determine the model of electrodiffusion in the junction. A 
physically sensible approximation to these reservoir solutions has been found 
to lead to radiation BCs for the Painleve II ODE which lies at the heart of 
the model, while the exact solutions lead to highly nonlinear BCs for that 
ODE. 

These new BCs, either exact or approximate, are not preserved in form 
under Backlund transformations, and it remains to classify BCs, as gener- 
ated by such transformations, that are physically relevant. The application 
of Backlund transformations to link boundary value problems of physical 
relevance does arise elsewhere, in particular in elasticity [12]. There mixed 
boundary value problems associated with normal loading and application of 
torsion may be linked in that way. Further applications involve boundary 
value problems describing indentation of shear-strained nonlinear elastic ma- 
terials [33]. We can be optimistic, therefore, that Backlund transformations 
can also be shown to relate physically interesting BCs in the present context. 

We have shown for the present model that Backlund transformations 
do leave invariant two quantities. One represents the total ionic (osmotic) 
pressure within the junction. The other represents the total resistive force per 
unit volume that the solute offers against ionic transport. We may hope that 
the invariance of these quantities provides a clue in the ongoing study of BCs 
and their relationship to flux-quantization by Backlund transformations. 



Appendix A: Gambier transformations of ex- 
act solutions 



Consider again the exact case of Planck's solution, as in (40). Note that for 
this solution, B as in ^ takes the value 2co, and A + = A_ = A, as required 
for application of G+- We consider a small perturbation of this solution in 
the form 



c±(x) = c + Ax + e d±(x) , E(x) = eF(x), A± = A. 



(Al) 



where d±(x) and F(x) are to be determined by requiring that c±(x), E(x) 
and A± satisfy ^ to first-order in e, given that they satisfy the system to 
zeroth-order because (40) is already a solution. Substituting from (Al) into 
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^ gives at first-order 

d + '(x) = (c + Ax) F(x) , gL '(x) = -(co + Ax) , 

X 2 F'(x) = d+(x) - d-(x) , 

and hence 

X 2 F"{x) = 2{c + Ax) F{x) . 
This has the general solution 

F{x) = oAi(s) + 6Bi(s) , s = 2(c + Ac)/(4A 2 A 2 ) 1/3 
where Ai and Bi are Airy functions of the first and second kind 



(A2) 



(A3) 



(A4) 
and a, 

b are arbitrary constants. Substituting from (Al) into (26) after noting that 
c + (x) — c_(a?) = e[d + (x) — d_(x)] 



= e\ 2 F'(x) 

= e(JA 4 ) 1/3 [aAi / (s) + 6Bi , (s), 
and then letting e — > 0, we obtain the new solution 

c+(x) = , c_(x) = 2A 2 F'(a;) 2 /i 71 (a;) 2 - 4Ax - 4c, 



(A5) 



= 2F'(x)/F(x) , i+ = , A. 



-4A, B = -4c n 



(A6) 



with F and F' as in (A4) and (A5). The treatment of Q_ , from the same 



starting point, is similar, and leads to the conjugate solution to (A6), as 



follows from (28). It is easily checked also that, just as the inverse Gambier 



transformations carry the seed solution (40 ) for the sequence of exact rational 



solutions of PII into the seed solution for the sequence of exact Airy solutions, 
so the Gambier transformations Q± carry the seed solutions for these Airy 



sequences, namely (A6) and its conjugate, back into the seed solution (40) 



for the sequence of rational solutions. 

Remark: These calculations reveal, for the first time, a connection be- 
tween the approximate solution [5] of PII in terms of Airy functions, and 
the known sequences of exact solutions of PII in terms of Airy functions [11]. 



The former is essentially E(x) as in (Al ), (A4), while the latter are generated 



from (A6) and its conjugate by Backhand transformations. 
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Appendix B: Exact boundary conditions and 
reservoir solutions 



The exact solution of (62), with the behaviour at x = — oo described there, 
is known to be 130, 311 



<p(x) = 2 In [(1 + Ae x/Xo ) / (l - Ae x/Xo )] , 



-1<A<1 



(Bl) 



where Aq is as in (64). From (60) and (|6T|) we then get 



c + (x) = c^oo) [(1 - Ae x ^) I (1 + Ae x ^)] 2 , 

c-(x) = c hoo) [(l + Ae x / X ») I (1 - Ae x ' x «)] 2 , 

E(x) = -4Ae x/x °/X (1 - A 2 e 2x/Xo ) . (B2) 
throughout the left-hand reservoir. The approximations are recovered when 



exp(x/A ) <C 1, so that (66) hold as before. 



From (B2) we have in particular that 

c + (0) = c ( _ oo) (l-A) 2 /(l + A) 2 , 

c_(0) = c ( _ oo) (l + A) 2 /(l-A) 2 , 

£(0) = -4A/A (l-A 2 ), (B3) 

from which we get 

c+(0)c_(0) = c ( „ 00 ) 2 , \ E(0) = y/2c+(0) - y/2c-(0) . (B4) 

These are the nonlinear exact BCs that now apply at the left-hand reservoir- 
junction interface when solving (|5]). Similarly, at the right-hand interface we 
get 



c+(l)c_(l) = c (+oo) 2 , Ai£(l) = <y/2c+(l)-<y/2c_(l). 



(B5) 



With a little effort it can be seen that the conditions (B4) and (B5) can also 
be written in the form 



c±(0) = c ( _ oo) + 1A 2 £(0) 2 ± ±A£(0) ,/8c ( _ oo) + A 2 £(0) 2 



c±(l) = c (+oo) + ±A 2 £(1) 2 ± |A£(1) a/8c (+oo) + A 2 £(l) 2 



(B6) 
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from which the approximate forms ( 66 ) are easily recovered 



The nonlinear BCs that apply for PII in the slab in this case can be seen 



from (B6); nothing is known about existence or uniqueness of its solutions. 



Once again it can be checked that BCs of the form (B6) are not preserved 
during flux-quantization by Backlund transformations. 
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